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ABSTRACT 

We use the 2dF Galaxy Redshift Survey to measure the dependence of the bj-band 
galaxy luminosity function on large-scale environment, defined by density contrast in 
spheres of radius 8/i~^Mpc, and on spectral type, determined from principal compo- 
nent analysis. We find that the galaxy populations at both extremes of density differ 
significantly from that at the mean density. The population in voids is dominated 
by late types and shows, relative to the mean, a deficit of galaxies that becomes in- 
creasingly pronounced at magnitudes brighter than Mbj — 5 logj^Q — 18.5. In con- 
trast, cluster regions have a relative excess of very bright early-type galaxies with 
Mb J — 5 log^o h'^ — 21. Differences in the mid to faint-end population between envi- 
ronments are significant: at Mb, — Slogj^g ^ = early and late-type cluster galaxies 
show comparable abundances, whereas in voids the late types dominate by almost 
an order of magnitude. We find that the luminosity functions measured in all den- 
sity environments, from voids to clusters, can be approximated by Schechter functions 
with parameters that vary smoothly with local density, but in a fashion which differs 
strikingly for early and late-type galaxies. These observed variations, combined with 
our finding that the faint-end slope of the overall luminosity function depends at most 
weakly on density environment, may prove to be a significant challenge for models of 
galaxy formation. 

Key words: galaxies: statistics, luminosity function — cosmology: large-scale struc- 
ture. 
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1 INTRODUCTION 

The galEixy luminosity function has played a central role in 
the development of modern observational and theoretical as- 
trophysics, and is a well established and fundamental tool 
for measuring the large-scale distribution of galaxies in the 
universe (Efstathiou, Ellis & Peterson 1988; Lovcday et al. 
1992; Marzke, Huchra & Geller 1994; Lin et al. 1996; Zucca 
et al. 1997; Ratcliffe et al. 1998; Norberg et al. 2002a; Blan- 
ton et al. 2003a;). The galaxy luminosity function of the 2dF 
Galaxy Redshift Survey (2dFGRS) has been characterised 
in several papers: Norberg et al. (2002a) consider the survey 
as a whole; Folkes et al. (1999) and Madgwick et al. (2002) 
split the galaxy population by spectral type; De Propris et 
al. (2003) measure the galaxy luminosity function of clus- 
ters in the 2dFGRS; Eke et al. (2004) estimate the galaxy 
luminosity function for groups of different mass. Such tar- 
geted studies are invaluable if one wishes to understand how 
galaxy properties are influenced by external factors such as 
local density environment (e.g. the differences between clus- 
ter and field galaxies). 

A natural extension of such work is to examine a wider 
range of galaxy environments and how specific galaxy prop- 
erties transform as one moves between them, from the very 
under-dense 'void' regions, to mean density regions, to the 
most over-dense 'cluster' regions. In order to 'connect the 
dots' between galaxy populations of different type and with 
different local density a more comprehensive analysis needs 
to be undertaken. Although progress has been made in this 
regard on both the observational front (Bromley et al. 1998; 
Christlein 2000; Hiitsi et al. 2002) and theoretical front (Pee- 
bles 2001; Mathis & White 2002; Benson et al. 2003; Mo 
et al. 2004), past galaxy redshift surveys have been severely 
limited in both their small galaxy rmmbers and small sur- 
vey volumes. Only with the recent emergence of large galaxy 
redshift surveys such as the 2dFGRS and also the Sloan Dig- 
ital Sky Survey (SDSS) can such a study be undertaken with 
any reasonable kind of precision (for the SDSS, see recent 
work by Hogg et al. 2003, Hoyle et al. 2003, and Kauffmann 
et al. 2004). 

In this paper we use the 2dFGRS galaxy catalogue to 
provide an extensive description of the luminosity distribu- 
tion of galaxies in the local universe for all density envi- 
ronments within the 2dFGRS survey volume. In addition, 
the extreme under-dense and over-dense regions of the sur- 
vey arc further dissected as a function of 2dFGRS galaxy 
spectral type, rj, which can approximately be cast as early 
and late-type galaxy populations (Madgwick et al. 2002, see 
Section 2). The void galaxy population is especially inter- 
esting as it is only with these very large survey samples 
and volumes possible to measure it with any degree of accu- 
racy. Questions have been raised (e.g., Peebles 2001) as to 
whether the standard ACDM cosmology correctly describes 
voids, most notably in relation to reionisation and the sig- 
nificance of the dwarf galaxy population in such under-dense 
regions. 

This paper is organised as follows. In Section 2 we pro- 
vide a brief description of the 2dFGRS and the way in which 
we measure the galaxy luminosity function from it. The lu- 
minosity function results are presented in Section 3, and 
then compared with past results in Section 4. We discuss 
the implications for models of galaxy formation in Section 5. 



Throughout we assume a ACDM cosmology with parame- 
ters Qm = 0.3, Qa = 0.7, and Ho = 100/i"^kms"^Mpc"^ 



2 METHOD 

2.1 The 2dFGRS survey 

We use the completed 2dFGRS as our starting point (Colless 
et al. 2003), giving a total of 221,414 high quality redshifts. 
The median depth of the full survey, to a nominal magni- 
tude limit of bj 19.45, is « « 0.11. We consider the two 
large contiguous survey regions, one in the south Galactic 
pole and one towards the north Galactic pole. To improve 
the accuracy of our measurement our attention is restricted 
to the parts of the survey with high redshift completeness 
(> 70%) and galaxies with apparent magrutude bj < 19.0, 
well within the above survey limit (see also Appendix C). 
Our conclusions remain unchanged for reasonable choices of 
both these restrictions. Full details of the 2dFGRS and the 
construction and use of the mask quantifying the complete- 
ness of the survey can be found in Colless et al. (2001, 2003) 
and Norberg et al. (2002a). 

Where possible, galaxy spectral types are determined 
using the principal component analysis (PCA) of Madgwick 
et al. (2002) and the classification quantified by a spectral 
parameter, rj. This allows us to divide the galaxy sample 
into two broad classes, conventionally called late and early 
types for brevity. The late types are those with 77 > —1.4 
that have active star formation and the early types are the 
more quiescent galaxy population with rj < —1.4. Approxi- 
mately 90% of the galEixy catalogue can be classified in this 
way. This division at 77 = —1.4 corresponds to an obvious 
dip in the rj distribution (Section 2.4; sec also Madgwick 
et al. 2002) and a similar feature in the bj — rr colour distri- 
bution, and therefore provides a fairly natural partition be- 
tween early and late types. When calculating each galEixy's 
absolute magnitude we apply the spectral type dependent 
k + e corrections of Norberg et al. (2002a); when no type 
can bo measured we use their mean k + e correction. In this 
way all galaxy magnitudes have been corrected to zero red- 
shift. 



2.2 Local density measurement 

The 2dFGRS galEixy catalogue is magnitude limited; it has a 
fixed apparent magnitude limit which corresponds to a faint 
absolute magnitude limit that becomes brighter at higher 
redshifts. Over any given range of redshift there is a cer- 
tain range of absolute magnitudes within which all galaxies 
can be seen by the survey and are thus included in the cat- 
alogue (apart from a modest incompleteness in obtaining 
the galaxies' redshifts). Selecting galaxies within these red- 
shift and absolute magnitude limits defines a volume-Umited 
sub-sample of galaxies from the magnitude-limited catalogue 
(see e.g. Norberg et al. 2001, 2002b; Croton et al. 2004); this 
sub-sample is complete over the specified redshift and abso- 
lute magnitude ranges. 

To estimate the local density for each galaxy we first 
need to establish a volume-limited density defining popu- 
lation (DDP) of galaxies. This population is used to flx 
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the density contours in the redshift space volume contain- 
ing the magnitude-limited galaxy catalogue. We restrict the 
magnitude-limited survey to the redshift range 0.05 < z < 
0.13, giving an effective sampling volume of approximately 
7x 10^/i~^Mpc^. Such a restriction guarantees that all galax- 
ies in the magnitude range — 19 > Mb, — 5 logj^Q h > —22 (i.e. 
effectively brighter than M* + 0.7) are volume limited, and 
allows us to use this sub-population as the DDP. The mean 
number density of DDP galaxies is 8.6 x W~'^h^Mpc~^ . In 
Appendix B we consider the effect of changing the magni- 
tude range of the DDP and find only a very small difference 
in our final results. 

The local density contrast for each magnitude-limited 
galaxy is determined by counting the number of DDP neigh- 
bours within an 8/i~^Mpc radius, Ng, and comparing this 
with the expected number, Ng, obtained by integrating un- 
der the published luminosity 2dFGRS function of Norberg 
et al. (2002a) over the same magnitude range that defines 
the DDP: 



Pa 



Ng-Ng 

No 



(1) 



i?=8h-lMpc 



In Appendix B we explore the effect of changing this smooth- 
ing scale from between 4/i~^Mpc to 12/i~^Mpc. We find that 
our conclusions remain unchanged, although, not surpris- 
ingly, smaller scale spheres tend to sample the under dense 
regions differently. Spheres of 8/i~^Mpc are found to be the 
best probe of both the under and over-dense regions of the 
survey. 

With the above restrictions, the magnitude-limited 
galaxy sample considered in our analysis contains a total 
of 81, 387 (51, 596) galaxies brighter than Mb, — 5 logj^g ^ ~ 
-17 (-19), with 30,354 (23,043) classified as early types 
and 42, 772 (23, 815) classified as late types. Approximately 
70% of all galaxies in this sample are sufficiently within the 
survey boundaries to be given a local density. Details of the 
different sub-samples binned by local density and type are 
given in Table 1. 



2.3 Measuring the luminosity function 

The luminosity function, giving the number density of galax- 
ies as a function of luminosity, is conveniently approximated 
by the Schechter function (Schechter 1976, see also Norberg 
et al. 2002a): 



= (t)*(L/L*)°' exp(-L/L*) d{L/L*) 



(2) 



dependent on three parameters: L* (or equivalently M*), 
providing a characteristic luminosity (magnitude) for the 
galaxy population; a, governing the faint-end slope of the 
luminosity function; and 0*, giving the overall normalisa- 
tion. Our method, which we describe below, will be to use 
the magnitude-limited catalogue binned by density and type 
to calculate the shape of each luminosity function, draw on 
restricted volume-limited sub-samples of each to fix the cor- 
rect luminosity function normalisation, then determine the 
maximum likelihood Schechter function parameters for each 
in order to quantify the changing behaviour between differ- 
ent environments. 

The luminosity function shape is determined in the 
standard way using the step-wise maximum likelihood 
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Figure 1. A comparison of the published 2dFGRS luminosity 
function (circles and dotted lines) to that calculated by our joint 
SWML (shape) /CiC (normalisation) method (squares and trian- 
gles) for select galaxy samples. Shown are the (top) full catalogue 
luminosity function (Norberg et al. 2002a) and cluster galaxy pop- 
ulation luminosity function (De Propris et al. 2003), and (bottom) 
the luminosity function for early and late-type galaxy sub-samples 
separately (Madgwick et al. 2002). 



method (SWML Efstathiou, EUis & Peterson 1988) and 
the STY estimator (Sandage, Tammann & Yahil 1979). See 
Norberg et al. (2002a) for a complete description of the 
application of these two techniques to the 2dFGRS. All 
STY fits are performed over the magnitude range — 17 > 
Mb, -51ogioft > -22. 

Such techniques fail to provide the luminosity function 
normalisation, however, and one needs to consider carefully 
how to do this when studying galaxy populations in dif- 
ferent density environments. To normalise each luminosity 
function we employ a new counts in cells (CiC) technique 
which directly calculates the number density of galaxies as a 
function of galaxy magnitude from the galaxy distribution. 
Briefiy, this is achieved by counting galaxies in restricted 
volume-limited sub-regions of the survey. We discuss our 
CiC method in more detail in Appendix A. As we show 
there, when galaxy numbers allow a good statistical mea- 
surement the luminosity function shape determined by the 
SWML and CiC methods agree very well. As the SWML 
estimator draws from the larger magnitude-limited survey 
rather than the smaller CiC volume-limited sub-samples, 
we choose the above two-step SWML/CiC approach rather 
than the CiC method alone to obtain the best results for 
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Figure 2. The distribution of the spectral parameter, ri, for the void, mean, and cluster galaxies used in our analysis (Table 1). The 
vertical dotted line at r; = — 1.4 divides the quiescent galaxy population (early types on left) from the more active star-forming galaxies 
(late types on right). From void to cluster environment, the dominant galaxy population changes smoothly from late-type to early-type. 



each luminosity function. Once the CiC luminosity function 
has been calculated for the same galaxy sample, the SWML 
luminosity function is then given the correct amplitude by 
requiring that the number density integrated between the 
magnitude range —19 > Mb, — 5 logj^Q h > —22 be the same 
as that for the CiC result. 



2.4 Comparison to previous 2dFGRS results 

In Fig. 1 we give a comparison of our measured luminosity 
functions for selected galaxy populations with the equivalent 
previously-published 2dFGRS results (see each reference for 
complete details). These include (top panel) the full survey 
volume (Norberg et al. 2002a) and cluster galaxy luminosity 
functions (De Propris et al. 2003), and (bottom panel) the 
luminosity functions derived for late and early-type galaxy 
populations separately (Madgwick et al. 2002). For all, the 
squares and triangle symbols show our hybrid SWML/CiC 
values while the circles and dotted lines give the correspond- 
ing published 2dFGRS luminosity function data points and 
best Schechter function estimates, respectively. The close 
match between each set of points confirms that our method 
is able to reproduce the published 2dFGRS luminosity shape 
and amplitude successfully. 

There are a few points to note. Firstly, the cluster lu- 
minosity function is not typically quoted with a value of 
</>* since the normalisation of the cluster galaxy luminos- 
ity distribution will vary from cluster to cluster (dependent 
on cluster richness). Because of this we plot the De Propris 
et al. cluster luminosity function using our </!>* value. 

Secondly, the Madgwick et al. early and late-type galaxy 
absolute magnitudes include no correction for galaxy evolu- 
tion, which, if included, would have the effect of dimming 
the galaxy population somewhat. We have checked the sig- 
nificance of omitting the evolution correction when deter- 
mining the galaxy absolute magnitudes and typically find 
only minimal differences in our results and no change to our 
conclusions. 

Thirdly, the STY Schechter function values we mea- 
sure tend to present a slightly 'flatter' faint-end slope than 
seen for the full survey: our all-type STY estimate returns 
a = —1.05 ± 0.02 (Table 1) whereas for the completed sur- 
vey (across the redshift range 0.02 < z < 0.25) the recov- 



ered a value is —1.18 ± 0.02 (Cole et al. in preparation). 
This difference is due primarily to three systematic causes: 
the minimum redshift cut required to define the DDP which 
results in a restricted absolute magnitude range over which 
we can measure galaxies; the non-perfect description of the 
galaxy luminosity function by a Schechter function together 
with the existing degeneracies in the M* — a plane; the 
sensitivity of the faint-end slope parametrisation to model 
dependent corrections for missed galaxies. For our results, 
these systematic effects do not hinder a comparison between 
sub-samples, but it is essential to take into account the dif- 
ferent cuts we imposed for any detailed comparison with 
other works. In Appendix C we discuss these degeneracies 
and correlations further. We test their influence by flxing 
each Q at the published field value when applying the STY 
estimator and find a typical variation of less than 0.2 mag- 
nitudes in Al* from the main results presented in Section 3. 
Such systematics do not change our conclusions. 

Lastly, the 2dFGRS photometric calibrations have im- 
proved since earlier luminosity function determinations, and 
thus the good match seen in Fig. 1 demonstrates that the 
new calibrations have not significantly altered the earlier 
results. 

In Fig. 2 we show the rj distribution for our void, mean, 
and cluster galaxy samples. The mean galaxy distribution 
is essentially identical to that shown in Fig. 4 of Madg- 
wick et al. (2002) for the full survey, demonstrating that 
the mean density regions contain a similar mix of galaxy 
types to that of the survey as a whole. For under-dense re- 
gions late types progressively dominate, while the converse 
is true in the over-dense regions. This behaviour can be un- 
derstood in terms of the density-morphology relation (e.g. 
Dressier 1980), and will be explored in more detail in the 
next section. 



3 RESULTS 

3.1 Luminosity functions 

The top panel of Fig. 3 shows the 2dFGRS galaxy lumi- 
nosity function estimated for the six logarithmically-spaced 
density bins and additional extreme void bin, Sg < —0.9, 
given in Table 1. The luminosity function varies smoothly 
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Figure 3. (top) The SWML luminosity functions for the 2dFGRS galaxy catalogue in regions of the survey of varying density contrast, 

Sfit from void to mean density to cluster. The best-fit Schechter function parameters for each arc given in Table 1 and the corresponding 
Schechter function curves arc over plotted here with dotted lines, (bottom) The void and cluster luminosity functions normalised to the 
mean luminosity function so as to highlight the relative differences in the shape of each distribution. The solid lines and bounding dotted 
lines show the appropriate Table 1 Schechter functions normalised to the mean Schechter function and 1<t uncertainty. 



as one moves between the extremes in environment. Each 
curve shows the characteristic shape of the Schechter func- 
tion, for which we show the STY fit across the entire range 
of points plotted with dotted lines. The Schechter parame- 
ters are given in Table 1, along with the number of galaxies 
considered in each density environment and the volume frac- 
tion they occupy. A number of points of interest regarding 



the variation of these parameters with local density will be 
discussed below. 

To examine the relative differences in the void and clus- 
ter galaxy populations with respect to the mean, in the bot- 
tom panel of Fig. 3 we plot the ratio of the void and clus- 
ter luminosity functions to the mean luminosity function. 
Also shown is the ratio of the corresponding Schechter func- 
tions to the mean Schechter function (solid lines) and la 



6 Croton et al. 



Table 1. Properties of our magnitude-limited galaxy samples, split by spectral type (all, early and late) and in seven density ranges 
(defined by Ss^^^ and is^^^, the density contrast in spheres of radius 8?i~^Mpc). The all-type sample is also split into an 'extreme void' 
sample. Nqal and /vol are the number of galaxies in each density bin and the volume fraction these galaxies occupy, respectively, /vol 
is given for all galaxy types only: car ly/latc- type density populations are just sub-divisions of the associated all-type sample. M* and a 
are the likelihood estimated Schechter function parameters, and <p* the associated normalisation. The integrated luminosity density, as 
defined by Eq. 3 with Lmin = 0, is given in the last column. All errors on the derived parameters reflect only the associated statistical 
uncertainty. 



Galaxy 


Galaxy 


°mm 




Ngal 


/vol 


M* 


a 




(pl) 


Type 


Sample 










- 5 logio 






10®hLgMpc~3 


all types: 


full volume 


-1.0 


oo 


81, 387 


1.0 


-19.65 ± 0.02 


-1.06 ± 0.02 


21.3 ± 0.6 


2.10 ± 0.08 




extreme void 


-1.0 


-0.90 


260 


0.09 


-18.26 ± 0.33 


-0.81 ± 0.50 


3.17 ±0.90 


0.08 ± 0.04 




void 


-1.0 


-0.76 


1, 157 


0.20 


-18.84 ± 0.16 


-1.06 ± 0.24 


3.16 ± 0.66 


0.16 ± 0.04 






-0.75 


-0.43 


3,331 


0.19 


-19.20 ± 0.10 


-0.93 ± 0.11 


5.99 ± 0.62 


0.36 ± 0.06 




mean 


-0.43 


0.32 


11, 877 


0.30 


-19.44 ± 0.06 


-0.94 ± 0.05 


11.3 ± 0.7 


0.86 ± 0.07 






0.32 


2.1 


21, 989 


0.24 


-19.64 ± 0.04 


-0.99 ± 0.04 


22.9 ± 1.0 


2.16 ± 0.13 






2.1 


6.0 


15, 656 


0.07 


-19.85 ± 0.06 


-1.09 ± 0.04 


49.0 ± 3.0 


5.96 ± 0.49 




cluster 


6.0 


oo 


3, 175 


0.01 


-20.08 ± 0.13 


-1.33 ± 0.11 


60.7 ± 13.2 


11.6 ± 3.4 


late type; 


full volume 


-1.0 


oo 


42, 772 




-19.30 ± 0.03 


-1.03 ± 0.03 


16.0 ± 0.6 


1.06 ± 0.06 




void 


-1.0 


-0.76 


855 




-18.78 ± 0.19 


-1.14 ± 0.24 


2.42 ± 0.66 


0.11 ± 0.04 






-0.75 


-0.43 


2, 249 




-19.07 ± 0.12 


-0.95 ± 0.14 


4.64 ± 0.68 


0.26 ± 0.06 




mean 


-0.43 


0.32 


7,261 




-19.24 ± 0.07 


-1.00 ± 0.07 


8.03 ± 0.61 


0.53 ± 0.06 






0.32 


2.1 


11, 921 




-19.36 ± 0.06 


-1.04 ± 0.05 


15.5 ± 1.0 


1.17 ± 0.11 






2.1 


6.0 


7,596 




-19.37 ± 0.07 


-1.03 ± 0.07 


36.3 ± 2.9 


2.73 ± 0.31 




cluster 


6.0 


oo 


1,316 




-19.34 ± 0.18 


-1.09 ± 0.20 


54.0 ± 12.2 


4.09 ± 1.31 


early type: 


full volume 


-1.0 


oc 


30, 354 




-19.65 ± 0.03 


-0.65 ± 0.03 


8.80 ± 0.22 


0.75 ± 0.03 




void 


-1.0 


-0.75 


220 




-18.62 ± 0.33 


-0.15 ± 0.53 


0.67 ± 0.10 


0.02 ± 0.01 






-0.75 


-0.43 


861 




-19.16 ± 0.14 


-0.43 ± 0.24 


1.62 ± 0.17 


0.88 ± 0.02 




mean 


-0.43 


0.32 


3,873 




-19.38 ± 0.08 


-0.39 ± 0.11 


4.13 ± 0.19 


0.27 ± 0.02 






0.32 


2.1 


8,809 




-19.59 ± 0.06 


-0.62 ± 0.06 


10.6 ± 0.4 


0.84 ± 0.06 






2.1 


6.0 


7, 163 




-19.89 ± 0.06 


-0.81 ± 0.06 


24.2 ± 1.6 


2.67 ± 0.23 




cluster 


6.0 


oo 


1,731 




-20.13 ± 0.18 


-1.12 ± 0.14 


37.1 ± 7.7 


6.00 ± 1.76 



uncertainty (dotted lines, where only the error in M* and a 
has been propagated). For a non-changing luminosity func- 
tion shape, this ratio is a flat line whose amplitude reflects 
the relative abundance of the samples considered. For two 
Schechter functions differing in both a and M*, the faint- 
end of the ratio is most sensitive to the differences in a and 
the bright end to the differences in M* . We note that the 
error regions on the Schechter function flts shown here do 
not include the uncertainty of the mean sample, as the cor- 
relation of its error with the other samples is unknown. This 
panel reveals significant shifts in abundances at the bright 
end: in voids there is an increasing deficit of bright galax- 
ies for magnitudes Mbj — 51og/i< — 18.5, while clusters 
exhibit an excess of very luminous galaxies at magnitudes 
Mbj -51ogft< -21. 

It is well-established that early and late-type galaxy 
populations have very different luminosity distributions 
(Fig. 1). In Fig. 4 we explore the density dependence of 
these populations. The upper panels show the luminosity 
functions and their Schechter function fits, as in Fig. 3, but 
for (left panel) early types and (right panel) late types sep- 
arately. In the corresponding lower panels we show the ratio 
of each extreme density population to the mean density lu- 
minosity function, following the same format as the bottom 
panel of Fig. 3. (We note that the mean luminosity func- 
tions for each type used in this figure are both very similar 
in shape to that shown in the bottom panel of Fig. 1). The 
best-fit Schechter parameters are given in Table 1. Again we 
see a smooth transition in the galaxy luminosity function 
as one moves through regions of different density contrast. 



The lower left panel of Fig. 4 shows a significant variation 
of the bright end early-type galaxy population with respect 
to the mean, while at the faint end the changes are more 
ambiguous, but with Schechter fits that suggest some evo- 
lution into the denser regions. Note that, although the faint 
end of our early-type cluster Schechter function is primar- 
ily constrained by the mid- luminosity galaxies in the sample, 
our maximum liklihood Schechter parameters are cjuite close 
to that found by De Propris et al. (2003) for a comparable 
galaxy population but measured approximately three mag- 
nitudes fainter. In contrast to the early types, in the lower 
right panel of Fig. 4 late-type galaxies show little change 
in relative population between the mean and cluster envi- 
ronments and a possible "tilt" favouring the faint-end for 
low-density environments. Due to deteriorating statistics we 
do not consider the type dependent extreme void luminosity 
function which was introduced in Fig. 3. 

The essence of our results is best appreciated when we 
directly compare the early and late-type galaxy distribu- 
tions, separately for the cluster and void regions of the sur- 
vey, as shown in Fig. 5. This figure reveals a striking con- 
trast: the void population is composed primarily of medium 
to faint luminosity late-type galaxies, while for the cluster 
population early types dominate down to all but the faintest 
magnitude considered. This is the central result of our study, 
and shows the crucial role of accurately determining the am- 
plitude of the luminosity function, since the shape alone does 
not necessarily determine the dominant population of a re- 
gion. 
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Figure 4. Comparing both the (top) absolute and (bottom) relative distributions of (left) early-type galaxies in different density 
environments, and (right) late-type galaxies in different density environments. In the bottom panels the luminosity functions have again 
been normalised to the mean (each to their respective type) as done previously in Fig. 3 (note that the shape of the mean for each type 
is very similar to that shown in Fig. 1). Here the solid lines and bounding dotted lines show the appropriate Table 1 Schechter functions 
normalised to the mean Schechter function and Icr uncertainty. 



3.2 Evolution with environment 

It is well known that the Schechter function parameters 
are highly correlated. In Fig. 6 we show the la (68% 2- 
parameter) and 3cr (99% 2-parameter) contours in the 
AI* — a plane for the early-type, late-type, and combined 
type cluster and void populations. For a given spectral type, 
all show a greater than 3a difference in the STY Schechter 
parameters between the void and cluster regions. Intermedi- 
ate density bins are omitted for clarity but follow a smooth 
progression with smaller error ellipses between the two ex- 
tremes shown. In Appendix C we explore in more detail the 
M* — a degeneracy and confirm that our results are robust. 

Our findings show that the galaxy luminosity function 
changes gradually with environment. We quantify this be- 
haviour in Fig. 7 by plotting the variation of M* and a as 
a function of density contrast, where points to the left of 



^8 = represent the under-dense to void regions in the sur- 
vey, and points to the right of this are measured in the over- 
dense to cluster regions. Late-type galaxies display a con- 
sistent luminosity function across all density environments, 
from sparse voids to dense clusters, with a weak dimming of 
M* in the under-dense regions, and an almost constant faint- 
end slope. In contrast, the luminosity distribution of early- 
type galaxies differs sharply between the extremes in en- 
vironment: M* brightens by approximately 1.5 magnitudes 
going from voids to clusters, while the faint-end slope moves 
from Q ~ —0.3 in under-dense regions to around a ~ —1.0 
in the densest parts of the survey. 

Finally, in Fig. 8 we plot the mean luminosity per 
galaxy, {pL)/{pg), obtained by integrating the luminosity 
function for each set of Schechter function parameters from 
Table 1: 
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Figure 5. A direct comparison of tlie early and late-type galaxy 
populations in the cluster environment (top two luminosity func- 
tions) and void regions of the survey (bottom two luminosity 
functions). The void population is composed almost exclusively 
of faint late-type galaxies, while in the clusters regions the galaxy 



population brighter than Af^j - 
inantly of early types. 



■Slogjg h = —19 consists prcdom- 



4>{L) dL , {pl) = 



4){L)L dL 



(3) 



where Lmin is a somewhat arbitrary observational cutoff 
chosen at Mbj — 5 logjQ h = —17. This is both the limit down 
to which we confidently measure our luminosity functions, 
and also the hmit beyond which the Schechter function no 
longer provides a good fit to the early-type luminosity func- 
tion of Madgwick et al. (2002). The final column of Table 1 
gives the total luminosity density in the various density con- 
trast environments, computed by integrating the Schechter 
function with no cutoff, to allow easy comparison with past 
and future analyses; the contribution to the calculated {pl) 
from luminosities below the observational cutoff is less than 
a few percent. We note that (p^) is directly related to the 
density contrast, ^g, by definition. It is interesting to see 
that the early-type galaxies in Fig. 8 are, on average, about 
a factor of two brighter per galaxy than the late types, even 
though the late types dominate in terms of both number and 
luminosity density. For all galaxy populations, the mean lu- 
minosity per galaxy shows a remarkable constancy across 
the full range of density environments. 



4 COMPARISON TO PREVIOUS WORK 

Historically, work on the dependence of the luminosity func- 
tion on large scale environment has been restricted primar- 
ily to comparisons between cluster and field galaxies, due 
to insufficient statistics to study voids. (Note that 'field' 
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-18.5 -19.0 -19.5 -20.0 -20.5 
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Figure 6. The Icr (68% 2-parameter) and 3ct (99% 2-parameter) 
contours of constant in the M* — a plane for the void and 
cluster STY estimates (in each case the void fit is on the left, 
corresponding to a fainter M*). Galaxy types are identified by 
the linestyle given in the legend. Even at the 3(T level significant 
differences in the void and cluster Schechter function parameters 
for each galaxy type can be seen. 



samples are usually flux-limited catalogues which cover all 
types of environments.) One of the aims of this work is to 
elucidate the properties of galaxies in void environments and 
understand the relationship between clusters and voids. In 
this section, we briefly summarise previous observations and 
compare them with the results presented in Section 3. 

We have already shown in Fig. 1 and Section 2 that our 
cluster and field results are equivalent to the published 2dF- 
GRS luminosity function results of Norberg et al. (2002a), 
Madgwick et al. (2002), and De Propris et al. (2003). The 
latter paper explains their cluster luminosity function by 
demonstrating that the field luminosity function can be ap- 
proximately transformed into the cluster luminosity function 
using a simple model where the cluster environment sup- 
presses star formation to produce a dominant bright, early- 
type population (see Section 4.4 of their paper for details). 
We expand upon such models in the next section. 

Bromley et al. (1998) considered 18,278 galaxies in the 
Las Campanas Redshift Survey (LCRS) as a function of 
spectral type and high and low local density. We confirm 
(e.g.. Fig. 7) their qualitative finding that for early-type 
galaxies the faint-end slope steepens with density whereas 
late-type objects show little or no significant trend. We can- 
not make a quantitative comparison to their result because 
they do not give the definition of their low density sample. 

Hutsi et al. (2003) use the Early Data Release of the 
Sloan Digital Sky Survey (SDSS) and the LCRS to consider 
the galaxy luminosity function as a function of density field, 
but in two-dimensional projection so their results are not 
directly comparable to ours. They find a faint-end slope of 
a ~ —1.1 in all environments and an increase in M* of 
roughly 0.3 magnitudes between the under and over-dense 
portions of their data. This is broadly consistent with the 
more detailed results obtained here with the full 2dFGRS 
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Figure 8. The mean luminosity per galaxy as a function of den- 
sity environment for each galaxy type, calculated from Ea.l3lusing 
the Schechter function parameters given in Table 1. The integral 



is performed by choosing L„ 



so that Mrr 



5 logi 
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(Section 3.2). The mean luminosity per galaxy of early types is 
consistently about a factor of two brighter than their late type 
counterparts across all density environments. 



Figure 7. The maximum liklihood Schechter function M* and a 
parameters for each of the six density contrast regions in Table 1 
(Figs. 3 and 4). Each panel shows the result for individual samples 
split by spectral type (early/late) and both types combined. 



catalogue when one averages over our two most underdense 
bins and three most over-dense. In separate work, these au- 
thors also consider the environmental dependence of cluster 
and supercluster properties in the SDSS and LCRS (Einasto 
et al. 2003a, 2003b). They show an almost order of magni- 
tude increase in the mean cluster luminosity between ex- 
tremes in density (defined in two dimensions by smoothing 
over a projected 10/i~^Mpc radius around each cluster). A 
comparison of their results to ours (i.e. Fig. 7) suggests a 
correlation between galaxy, galaxy group, and galaxy cluster 
properties in a given density environment. A more detailed 
exploration would shed light on the connection between viri- 
alised objects of different mass with local density. We defer 
such an investigation to later work. 

In a series of papers, members of the SDSS team under- 
took an analysis of the properties of galaxy samples drawn 
from under-dense regions in the SDSS (Rojas et al. 2004, 
2003; Goldberg et al. 2004; Hoyle et al. 2003). Of most rele- 
vance to our study is the work of Hoyle et al. who completed 
a preliminary analysis of the SDSS void luminosity function, 
defined in regions of Sj < —0.6 using a smoothing scale of 
7/i~^Mpc. Their sample of 1, 010 void galaxies are typically 
fainter and bluer than galaxies in higher density environ- 
ments but with a similar faint-end slope. Their results are 



consistent with what we find using a sample which contains 
about twice the number of void galaxies as defined by Hoyle 
et al.. Using the same void galaxy catalogue, Rojas et al. 
(2003) show that this behaviour is not merely an extrapo- 
lation of the density-morphology relationship (e.g. Dressier 
1980) into sparser environments. By measuring the concen- 
tration and Sersic indices (Sersic 1968, Blanton et al. 2003b) 
of void and field galaxies they detect no significant shift in 
the morphological mix, even though their void galaxy sam- 
ple is shown to be significantly bluer. 

Also using the SDSS dataset, Hogg et al. (2003) consider 
the mean environment as a function of luminosity and colour 
of 115, 000 galajdes, on smoothing scales of 1 and 8/i~^Mpc. 
They find that their reddest galaxies show strong correla- 
tions of luminosity with local density at both the faint and 
bright extremes, whereas the luminosities of blue galaxies 
have little dependence on environment. These conclusions 
are consistent with the present results for our early (red) and 
late (blue) type luminosity functions (Fig. 7). However by 
restricting attention to the average environment of a galaxy 
of given luminosity and colour, their sample is by defini- 
tion dominated by galaxies in over-dense environments. The 
measures they consider are therefore insensitive to one of the 
main questions of interest to us here, namely whether the 
characteristic galaxy population in the voids is distinctively 
different from that in other density environments. Indeed, 
we clearly find evidence for a population which is particu- 
larly favoured in void regions, namely faint late- type galaxies 
(Fig. 5). 
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Table 2. A mnemonic summary of our main results, drawing on the work of De Propris et al. (2003) and Mo et al. (2004) to interpret 
the observed behaviour in Figs. 5 and 7 in terms of physical processes which govern the void and cluster galaxy populations. 



Region 


Observation 


Process 


Voids 


faint, late type 
galaxies dominate 


1. galaxies typically reside at the centre 
of low mass dark halos {=> faint) 

2. gas is available for star formation (=i> blue) 

3. merger rate is low (^ spirals) 


Clusters 


mid-bright, early-type 
galaxies dominate 


1. typically satellite and central galaxies of 
massive dark halos (=^ mid-bright) 

2. gas is unavailable for star formation (=^ red) 

3. merger rate is high (=> ellipticals) 



5 DISCUSSION 

As clusters are comparatively well-studied objects, and have 
already been addressed using the 2dFGRS by De Propris et 
al. (2003), we focus here primarily on a discussion of the 
voids. 

A detailed analysis of void population properties has 
recently become possible due to significant improvements 
in the quality of both theoretical modelling and observa- 
tional data, as summarised by Benson et al. (2003). Pee- 
bles (2001) has argued that, visually, observed voids do not 
match simulated ones and discussed several statistical mea- 
sures for quantifying a comparison, primarily the distance to 
the nearest neighbour in a reference sample. However the cu- 
mulative distributions of nearest neighbour distances shown 
in Figs. 4-6 of Peebles (2001) show very little difference be- 
tween the reference-reference and test-reference distribu- 
tions. It is not surprising that these statistical measures are 
insensitive to a void effect, since they are dominated by clus- 
ter galaxies. Our method is designed to overcome this diffi- 
culty by explicitly isolating the void population of galaxies 
so that their properties can be studied. 

Motivated by the claims of discrepancies in Peebles 
(2001), Mathis & White (2002) investigated the nature of 
void galaxies using N-body simulations with semi-analytic 
recipes for galaxy formation. They call into question the as- 
sertion of Peebles (2001) that ACDM predicts a population 
of small haloes in the voids, concluding that "the popula- 
tion of faint galaxies... does not constitute a void popula- 
tion". More specifically, they find that all types of galaxies 
tend to avoid the void regions of their simulation, down to 
their resolution limit of Mb = —16.27 in luminosity and 
Mb = —18.46 in morphology. 

The abundance of faint galaxies we find in the void re- 
gions of the 2dFGRS seems to be at odds with the Mathis & 
White predictions. However their results rely on the Peebles 
(2001) cumulative distribution of galaxies as a function of 
over density (their Fig. 3) which, like cumulative distribu- 
tions in general, are rather insensitive to numerically-minor 
components of the galaxy population. Note that Mathis & 
White define density contrast using the dark matter mass 
distribution smoothed over a 5 h^^Mpc sphere, whereas we 
measure the density contrast by galaxy counts. Another pos- 
sible source of discrepancy is the uncertainties in their semi- 
analytic recipes, such as the implementation of supernova 
feedback, which can strongly effect the faint-end luminosity 
distribution. 



There has been recent discussion in the literature about 
the nature of the faint-end galaxy population and its depen- 
dence on group and cluster richness. Most notably, TuUy et 
al. (2002) show a significant steepening in the faint-end pop- 
ulation as one considers nearby galaxy groups of increasing 
richness, from the Local Group to Coma. On the surface of 
it, this might seem at variance with our results, which are 
rather better described by a faint-end slope which is approx- 
imately constant with changing density environment for the 
full 2dFGRS galaxy sample (Fig. 7). However the steepening 
of the faint-end slope they find primarily occurs at magni- 
tudes fainter than Mb = —17, which is beyond the limit 
we can study with our sample. Also, their analysis focuses 
on individual groups of galaxies, while we have chosen to 
work with a much bigger galaxy sample and have smoothed 
it over a scale much larger than the typical cluster. Indeed, 
as discussed in Appendix B, when the smoothing scale is 
significantly larger than the characteristic size of the struc- 
tures being probed it is possible that the Schechter function 
parameters may become insensitive to the small scale shifts 
in population. This effect, of course, would be less significant 
for survey regions which host clusters of clusters (i.e. super- 
clusters), and which axe prominently seen in the 2dFGRS 
(Baugh et al. 2004, Croton et al. 2004). When sampling the 
2dFGRS volume the trend with density that one sees using 
4h^^ Mpc spheres in Fig. Bl is consistent with the TuUy 
et al. result, although one needs to additionally understand 
the influences of Poisson noise. 

TuUy et al. attribute their results to a process of pho- 
toionisation of the IGM which suppresses dwarf galaxy for- 
mation. Over-dense regions, which at later times become 
massive clusters, typically collapse early and thus have time 
to form a dwarf galaxy population before the epoch of reion- 
isation. Under-dense regions, on the other hand, begin their 
collapse at much later times and are thus subject to the 
photoionisation suppression of cooling baryons. This, they 
argue, explains the significant increase between the dwarf 
populations of the Local Group (low density environment) 
and Coma (over-dense environment). Although suggestive, 
a deeper understanding of what is happening will require a 
much more statistically significant sample. 

Recently Mo et al. (2004) have considered the depen- 
dence of the galaxy luminosity function on large scale envi- 
ronment in their halo occupation model. In this model, the 
mass of a dark matter halo alone determines the proper- 
ties of the galaxies. They create mock catalogues built with 
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a halo-conditional luminosity function (Yang et al. 2003) 
which is constrained to reproduce the overall 2dFGRS lu- 
minosity function and correlation length for both luminos- 
ity and type. They analyse their data by smoothing over 
spheres of radius 8 /i~^Mpc in their mock catalogue, and 
measure the luminosity function as a function of density 
contrast. Their work is performed in real space while we 
are restricted to work in redshift space. Nonetheless, their 
predictions qualitatively match our density-dependent lu- 
minosity functions; a quantitative comparison is deferred to 
subsequent work. 

In the framework of the Mo et al. model, the reason 
the faint-end slope a has such a strong dependence on lo- 
cal density for early types (Fig. 7) is that faint ellipticals 
tend to reside predominantly in cluster-sized halos. The a 
dependence is weaker for late types because faint later-type 
galaxies tend to live primarily in less massive haloes, which 
are present in all density environments. The correlations be- 
tween dark halo mass and the properties of the associated 
galaxies are not a fundamental prediction of their model, 
but are input through phenomenological functions adjusted 
to give agreement with the 2dFGRS overall luminosity func- 
tions by galaxy type. However it would be a non-trivial result 
if the correlations are the same independently of whether the 
dark matter halo is in a void or in a cluster. For instance, 
this property would not apply in models for which reionisa- 
tion more efficiently prevented star formation in under-dense 
environments than in over-dense environments, as discussed 
above. 

An interesting consequence of the Mo et al. (2004) 
model is that the luminous galaxy distribution (which is easy 
to observe but hard to model) correlates well with the dark 
halo mass distribution (which is hard to observe but easy to 
model) . If their predictions prove to give a good description 
of the present data it will lend credence to the underlying 
assumption of their model — that the environmental depen- 
dence of many fundamental galaxy properties are entirely 
due to the dependence of the dark halo mass function on 
environment. Exactly why this is would still need to be ex- 
plained, however such a demonstration may facilitate more 
detailed comparisons between theory and observation than 
previously possible. 

An important result of our work is presented in Fig. 5, 
where a significant shift in the dominant population between 
voids and clusters is seen. Such a result points to substan- 
tial differences in the evolutionary tracks of cluster and void 
early-type galaxies. Cluster galaxies have been historically 
well studied: they are more numerous and much brighter 
on average, with an evolution dominated by galaxy-galaxy 
interactions and mergers. In voids, however, the picture is 
not so clear. A reasonable expectation would be that the 
dynamical evolution of void galaxies should be much slower 
due to their relative isolation, with passively evolved stellar 
populations and morphologies similar to that obtained dur- 
ing their formation. Targeted observational studies of void 
early-type galaxies may reveal much about the high redshift 
formation processes that go into making such rare objects. 

Table 2 summarises our main results and provides a 
qualitative or mnemonic interpretation based on our obser- 
vations and the work of Mo et al. (2004) and De Propris et al. 
(2003) (and references therein). Our primary result is the 
striking change in population types between voids and clus- 



ters shown in Fig. 5: faint, late-type galaxies are overwhelm- 
ingly the dominant galaxy population in voids, completely 
the contrary of the situation in clusters. The existence of 
such a population in the voids, and more generally the way 
populations of difi'erent type are seen to change between dif- 
ferent density environments, place important constraints on 
current and future models of galaxy formation. 
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APPENDIX A: THE COUNTS IN CELLS 
LUMINOSITY FUNCTION ESTIMATOR AND 
COMPARISON TO THE SWML RESULTS 

Our counts-in-cells (CiC) method to measure the density 
dependent luminosity function and obtain its amplitude is 
simple and will be illustrated with the example of a mock 
galaxy sample in a cubical volume of side-length L. The full 
luminosity function for such a sample is trivial. By defini- 
tion it is simply the number of galaxies in each magnitude 
interval divided by the volume of the box; 

$(M) = N{M) I . (Al) 

To determine the luminosity function as a function of 
local galaxy density we require two additional pieces of in- 
formation. Firstly we sub-divide the galaxy population into 
density bins. The local density for each galaxy is calculated 
within an 8/i~^Mpc radius as described in Section 2.2. This 
gives us the number of galaxies in each density bin belonging 
to each magnitude range, Ns^{M). 

Secondly we determine the volume which should be at- 
tributed to the various density bins. We do this by finding 
the fraction of the volume in which the galaxies of each 
density bin reside, /^^ . This fraction is measured by mas- 
sively oversampling the box with randomly placed 8/i~^Mpc 
spheres, in each of which we estimate a local density in the 
same way as before. Once all spheres have been placed we 
count the number which have a local density in each density 
range. The volume fraction of each bin is then just the frac- 
tion of spheres found in each bin. Since the total volume of 
the box is known, the volume of each density bin is now also 
known. The density dependent luminosity function is then 
calculated as: 




-17 -18 -19 -20 -21 -22 
^^bj - Slog, oh 



Figure Al. A comparison of the raw counts in cells luminosity 
function with the normalised SWML luminosity function, as de- 
scribed in the text (Section 2.3 and Appendix A). Shown are the 
cluster, mean, and void populations consisting of all galaxy types 
only, although all luminosity functions used in this paper behave 
equally as well. The shapes estimated by the two very different 
methods are in very good agreement over the magnitude ranges 
considered. 



<l,s,{M) = Ns,{M) / fs,L\ (A2) 

The situation becomes more complicated when dealing 
with a magnitude limited redshift survey instead of a sim- 
ple simulated box. Galaxy counting and volume estimation 
must now be restricted to regions of the survey in which the 
magnitude range being considered is volume limited. This 
range of course changes for each magnitude bin in which 
the luminosity function is measured. In addition, small cor- 
rections (< 10%) are required when counting galaxies to 
account for the spectroscopic incompleteness of the survey 
(see Croton et al. 2004). In all other respects, however, the 
calculation of (M) is the same as in the "box" example 
given above. 

In Fig. Al we show a comparison of the 2dFGRS CiC 
and SWML luminosity functions calculated from the same 
void, mean, and cluster galaxy samples. The SWML lumi- 
nosity function has been normalised to the CiC measurement 
as described in Section 2.2. We see that both methods pro- 
duce almost identical luminosity distribution shapes. This 
gives us confidence that the CiC luminosity function can 
be used to normalise the SWML luminosity function in an 
unbiased way. 

Because of the volume-limited restriction of the CiC 
method, the number of galaxies used to calculate the lu- 
minosity function is smaller than for the SWML method, 
which draws from the larger magnitude-limited catalogue. 
However the benefit of the CiC method is that it gives a di- 
rect measurement of the number density of galaxies rather 
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Figure Bl. The difference in ttie STY Schecliter function pa- 
rameters when the local density is calculated with an increasingly 
fainter density defining population (DDP) : Mj^in — 5 log^Q h = 
—20, —19, and —18. Such a change also changes the redshift range 
of the included volume as described in the text. For clarity only 
results for all galaxy types are shown. The reference sample is the 
^min — Slogj^o ft- = — 19 DDP used throughout this paper, and 
the other DDP results are shown relative to this. 



Figure B2. The difference in the STY Schechter function pa- 
rameters for different density bins when calculated by smooth- 
ing the local galaxy distribution on different scales: 4, 8, and 
12ft.~^Mpc. For clarity only results for all galaxy types are shown. 
The reference sample is the 8/i~^Mpc sphere smoothing scale used 
throughout this paper, and the other smoothing scale results are 
shown relative to this. 



than just the shape of the distribution as the SWML estima- 
tor docs. In addition, the CiC method is very easy to apply 
to mock catalogues, as described above. By combining the 
CiC and SWML methods we capture the best features of 
both. 



APPENDIX B: THE EFFECTS OF CHANGING 
THE DENSITY DEFINING POPULATION AND 
SMOOTHING SCALE 

In our analysis we are required to make two important 
choices before beginning. The first is to find the widest pos- 
sible absolute magnitude range for the density defining pop- 
ulation (DDP, see Section 2.2) while maximising the amount 
of the 2dFGRS survey volume sampled. The second is the 
scale over which wo smooth the DDP galaxy distribution to 
determine the density contours within this volume. We will 
now consider the effect of changing each of these choices in 
turn. 

The density defining population 

The DDP is important in that it not only sets the mean 
density of galaxies used to define the density contours, but 
also determines the redshift range of the full magnitude- 
limited catalogue to be included in the analysis. Clearly 
one would like as high-statistics a sample as possible in as 
large a volume as possible for the best results. In a vol- 
ume limited galaxy sample such as the DDP, the maxi- 
mum galaxy redshift available is constrained by the spec- 
ified faint absolute magnitude limit: galaxies beyond this 



redshift range arc no longer guaranteed to be volume lim- 
ited and are thus not included. For the DDP faint magni- 
tude limit of Mmin — 5 logj^o h = —19 the maximum survey 
boundary is z = 0.13. Changing the faint magnitude limit to 
Mmin — 5 logj^Q /i = — 18 (20), i.e. a denser (sparser) DDP, re- 
sults in a maximum redshift boundaxy oi z = 0.088 (0.188), 
i.e. a smaller (larger) sampling volume. 

In Fig. Bl we show the result found when repeating 
the analysis of Section 3 (Fig. 7) but using DDPs defined 
by different faint absolute magnitude limits. We plot the 
STY M* and a values for each density bin relative to the 
Mmin — Slogio /i = — 19 DDP used throughout this paper. 
The faintest DDP shown, Mmin — 5 logj^Q h = —18, is approx- 
imately 8 times denser than the brightest, Mmi„ — 5 log^Q h = 
—20, but with a volume roughly 30 times smaller. Even so, 
almost all measurements shown across all density bins are 
consistent at the la level, demonstrating that our defini- 
tion of the DDP is a robust representation of the underlying 
global density distribution. 



The smoothing scale 

Now let us examine how changing the smoothing scale with 
which we define local density affects the shape of our lumi- 
nosity functions. In Fig. B2 we examine the values of the 
Schechter parameters when measured with spheres of radius 
4 and 12/i~'^Mpc, compared to when the luminosity function 
is measured with an 8/i~^Mpc sphere. 

Fig. B2 shows a typical deviation of < 0.2 magnitudes 
for M* and < 0.2 for a. The 4/i~'^Mpc smoothing scale de- 
viates strongly from the other values in the under-dense re- 
gions (the first two points lie beyond the axis range plotted) , 
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however in these environments such a smoothing scale gives 
a poor estimate of the local galaxy density due to Poisson 
noise in small number counts. Indeed, Hoyle et al. (2004) 
have shown that in the extreme under-dense 2dFGRS sur- 
vey regions the characteristic scale of voids is approximately 
15/i~^Mpc. For cluster regions 4/i~^Mpc spheres can be em- 
ployed and would give a higher resolution discrimination 
of the structure. The 8/i~^Mpc smoothing scale we have 
adopted captures the essential aspects of voids while roughly 
optimising the statistical signal, and is thus a good probe of 
both the under and over-dense regions of the survey volume. 



APPENDIX C: SYSTEMATIC EFFECTS WHEN 
ESTIMATING THE SCHECHTER FUNCTION 
PARAMETERS 

One may ask to what degree the trends seen in Fig. 7 and 
Table 1 are influenced by the systematics discussed in Sec- 
tion 2.4. We note there that our STY measurements recover 
a flatter faint end than the current published 2dFGRS lu- 
minosity function for the completed catalogue. We identify 
three systematic effects which contribute to this behaviour: 
the absolute magnitude range considered when applying the 
STY estimator, the fact that the luminosity function is not 
perfectly described by a Schechter function, and the sen- 
sitivity of the faint-end slope parameterisation to model- 
dependent corrections included to account for missed galax- 
ies. 

We find that the first two of these effects have the 
strongest influence on the measured STY faint-end slope. 
Indeed, testing the first reveals that any STY estimate of 
the 2dFGRS luminosity function over a restricted absolute 
magnitude range displays a systematic shift in the recov- 
ered STY parameters along a line in the M* — oi plane. 
The brighter the faint magnitude restriction, the flatter the 
faint-end slope is and the fainter the characteristic mag- 
nitude becomes. Such behaviour is a consequence of small 
but important deviations in the galaxy luminosity function 
shape from the pure Schechter function assumed by the STY 
estimator. In addition. Fig. 11 of Norberg et al. (2002) re- 
veals a dip in the luminosity function between the magni- 
tude range —17 > Mbj — 51ogjQ h > —18 and a steepening 
faintward of this. In our analysis only galaxies brighter than 
A/bj — 5 logj^Q h = —17 are considered due to the restriction of 
the DDP. This limitation adds extra weight to the influence 
of the dip on the STY fit contributing further to a fiatter 
estimation of a. When mock galaxy catalogues constructed 
to have a perfect Schechter function luminosity distribution 
are analysed in an identical way to the 2dFGRS samples, 
we find that the above systematics all but disappear and 
the "true" M* and a values are recovered for any reason- 
able choice of STY fitting range. 

Finally, we note that the sensitivity of the faint-end 
slope parametrisation to systematic corrections for spectro- 
scopically missed galaxies is minimised by restricting our 
analysis to galaxies with bj < 19, for which the spectro- 
scopic incompleteness is typically less than ~ 8% (see Fig. 16 
of CoUess et al. 2001). 

Given that a full correction of the above systematic ef- 
fects is not possible in our analysis, the next best thing we 
can do is try to quantify to what degree they influence our 
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Figure CI. The shift in the STY Schechter function parameter 
M* when a is kept fixed at the published field value, compared 
with that found when a remains free (Table 1 and Fig. 7). The 
points are plotted as a function of local density and shown for 
each galaxy type and the combined all-type sample. 



results and conclusions. We do this by fixing the faint-end 
slope a when applying the STY estimator: at —1.2 for the 
all- types samples, —1.1 for the late- type samples, and at 
—0.5 for the early-type samples. Such choices enforce the 
published field luminosity function faint end values found 
by Norberg et al. (2002a) and Madgwick et al. (2002) and 
remove the degeneracy in the M* — a plane. 

Fig. CI shows the size of the shift in M* when such 
constraints are applied relative to that found when a is al- 
lowed to remain free (i.e. Table 1 and Fig. 7). Most notable 
here is that, apart from the two most over-dense bins in the 
early-type sample, there is no significant difference in the be- 
haviour oi M* with density environment. The approximate 
0.2 magnitude offset seen in this figure can be understood by 
remembering that because the faint-end slope we measure 
when a remains free is slightly fiatter than the published val- 
ues (due to the systematics discussed above), by artifically 
fixing a one forces M* to move to compensate. The impor- 
tant point is that the trends seen in Fig. 7 with changing 
local density remain unchanged. 

For the two most over-dense early-type samples a shift 
of up to 0.4 magnitudes is seen. We note from Table 1 
that our best-fit (free a) early-type cluster value of q = 
-1.12 ± 0.14 is weU matched by the equivalent 2dFGRS De 
Propris et al. (2003) resuh of -1.05±0.04. In effect, by con- 
straining the early-type cluster faint-end slope to the field 
value of —0.5 we ignore the real changes in galaxy popula- 
tion seen between the Madgwick et al. (their Fig. 10) and 
De Propris et al. (their Fig. 3) luminosity functions (see also 
our Fig. 2). Such population changes, we argue, result in the 
strikingly different Schechter function parameterisation be- 
haviour seen in Fig. 7 and Table 1 for early and late-type 
galaxies. Fig. CI gives us confidence that the M* — a de- 
generacies and systematics investigated here are not signif- 
icantly infiuencing our results or the conclusions we draw 
from them. 



